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Abstract. A theory of higher rank multircsolution analysis is 
given in the setting of abclian multiscalings. This theory enables 
the construction, from a higher rank MRA, of finite wavelet sets 
whose multidilations have translates forming an orthonormal ba- 
sis in L 2 (R d ). While tensor products of uniscaled MRAs pro- 
vide simple examples we construct many nonscparablc higher rank 
wavelets. In particular we construct Latin square wavelets as rank 
2 variants of Haar wavelets. Also we construct nonseparable scal- 
ing functions for rank 2 variants of Meyer wavelet scaling functions, 
and we construct the associated nonseparable wavelets with com- 
pactly supported Fourier transforms. On the other hand we show 
that compactly supported scaling functions for biscaled MRAs are 
necessarily separable. 



1. Introduction 

The term multiscaling in wavelet theory commonly refers to the var- 
ious scaling levels present in a nest of subspaces 

• • • C V-! c v c v x c . . . 

of a multiresolution analysis in L 2 (M. d ). See, for example, Dutkay and 
Jorgensen [5]. However such subspaces are associated with the powers 
of a single dilation matrix. In contrast we develop here a theory of 
wavelets for higher rank multiresolution analyses which are generated 
by several independent commuting dilation matrices. 

Recall that a wavelet set, or multiwavelet, is generally taken to be a 
set of functions ipi(x), . . . , ipt{x) in L 2 (M. d ) for which certain translates 
of dilates form an orthonormal basis of the form 

{(det A)-^ipi(A m x + k) : m e Z, k e Z d , 1 < % < t}, 

where A is a scaling matrix in GL(M. n ). Wavelet theory is concerned 
with identifying constructions for which the wavelets ipi exhibit forms 
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of directionality and smoothness. In particular it has been of interest 
to obtain multivariable wavelets which are in some sense nonseparable 
with respect to the coordinates of R d . See for example [T], [3], [7], |9]. 

, | I12| , | |16] . In particular Belogay and Wang [1] construct nonsep- 
arable wavelets in M 2 , for some dilation matrices with determinant 2, 
which are arbitrarily smooth. 

Most commonly, particularly in multiresolution analysis, the dilation 
group is singly generated, as above. Some recent studies with multi- 
scalings that go beyond this are summarised in Gu, Labate, Lim, Weiss 
and Wilson |5|. In these settings volume-preserving sheering unitaries 
combine with a single strict dilation to generate the nonabelian dila- 
tion groups of interest, and associated wavelet sets are constructed. 
Also these so-called affine orthonormal systems exhibit specific forms 
of directionality and nonseparability. In contrast to this our wavelets 
are constructed for a dilation representation of the abelian group 27 '. 

There are two main approaches to constructing wavelet sets, namely 
the multiresolution analysis approach of Mallat |13| , on the one hand, 
and constructions associated with self-similar tilings of M. d on the other. 
See Wang | |16] , for example, for connections with tilings. As far as the 
authors are aware there has been little development of traditional mul- 
tiresolution analysis wavelet theory for multiscaled settings in which 
the role of dilation matrices A m is played by an abelian group of dila- 
tion matrices A™ 1 A™ 2 . . . A™ s . The simplest such multiscaled context 
of this kind is the biscaled dyadic case for wavelets in L 2 (IR 2 ) associated 
with the dilation matrices 



A = 



'2 


0" 


,B = 


"1 


0" 
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If iPa( x ) an d iPb{z) are univariate dyadic wavelets in L 2 (IR) then the 
separable function i/j(x, y) = ipA(x)ipB(y) is a biscaled wavelet, that is 
the set 

{2- (rn+n)/2 ij(A m B n x + k) : (m, n) e Z 2 , k e Z 2 }, 

is an orthonormal basis in L 2 (R 2 ). 

The evident nonlocality of such separable wavelets has possibly not 
encouraged the elaboration of a multiscaled wavelet theory. (See, for 
example, the discussion in |17j .) However, we shall show that even in 
this apparently adverse setting of separated coordinates it is possible to 
construct nonseparable wavelets and even nonseparable scaling func- 
tions. Moreover we develop what might be termed a theory of higher 
rank wavelets and multiresolution analysis. 
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We define a higher rank multiresolution analysis (0,21), of rank r, 
where 23 = {Vi : i G Z r } is a commuting lattice of closed subspaces of 
L 2 (M. d ) with appropriate inclusions. See Definition 12.31 and Section 2.8. 
In our first main result Theorem 13.41 we show how one may construct 
multiscaled wavelet sets from a higher rank multiresolution of rank two. 

The distinguished novelties of higher rank multiresolution analysis 
are already present in the simplest setting of biscaling (r = 2) and 
dimension 2 (d = 2). For such a biscaled multiresolution analysis 
(BMRA) the scaling function (p(x, y) possesses two marginal filter func- 
tions JTif (£)> m df(0 an d these must satisfy the intertwining relation 

(1) rni(BZ)m»(Z)=m$(Z)m*{A£)- 
This filter relation follows from the coincidence of the tripe subspace 
inclusions associated with the dilation pairs B, AB and A, BA. That 
is, the intertwining relation is a consequence of the lattice structure 
of 03. In addition to this the orthogonality structure, or commuting 
(projection lattice) structure, of the BMRA leads to a filter identity of 
some complexity namely, 



™£(£i + 7T, 2£ 2 )m£(£i, £ 2 )m?(£i, 6)™?(2£i, £ 2 + tt)- 
m£(&, 2&)m£(6 + 7T, 6X(6 + vr,6X(26,6 + 7r)- 
(£i + vr, 2£ 2 )m£(£ 1 , £ 2 + vr)mf (6, 6 + 7r)mJ(26, 6)+ 
m£(£i, 2£ 2 )m£(£i + 7T, 6 + 7r)m|(6 + 7r,6 + 7r)m|(26,6) = 0. 

These two necessary conditions present challenges for the construc- 
tion of nonseparable scaling functions. Nevertheless, in another of our 
main results, Theorem 17. 2\ we construct examples of BMRAs for which 
the scaling function is indeed nonseparable. In this case the derived 
wavelets are similarly nonseparable. Furthermore, they inherit smooth- 
ness from <fi. In this way we obtain singleton bidyadic wavelets whose 
Fourier transforms have compact support and these wavelets are in fact 
higher ranks variants of the well-known Meyer wavelets |14j . 

On the other hand, even for separable rank 2 multiresolutions and 
separable scaling functions we show that there is sufficient freedom in 
the nondyadic case 

(detA - l)(det£- 1) > 2 

to derive nonseparable wavelets. This is achieved by constructing filter 
matrix functions which are not elementary tensor products. The con- 
struction of these filter functions parallels the well-known method of 
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unitary matrix completion although our arguments, given in the proof 
of Theorem 13.41 involve a nesting of several Gram Schmidt completion 
processes. The wavelet sets here include some very interesting and 
computable examples of what we term Latin square wavelet sets, for 
evident reasons. (See Theorem 14.31 ) These are multiscaled versions of 
the classical wavelets associated with Haar bases. 

In another of our main results we reveal a striking constraint for 
compactly supported scaling functions in the biscaled theory in L 2 (IR 2 ), 
namely that such functions are necessarily separable. Equivalently put, 
a BMRA (0, 53) with compactly supported scaling function is equiv- 
alent to an elementary tensor of two uniscaled MRAs. This fact may 
have been a further implicit obstacle to the development of multiscaled 
multiresolution wavelets. 

There are many intriguing wavelet directions that now seem to beckon. 
For example, we have not particularly addressed smoothness and ap- 
proximation properties in this article. Although we have shown that 
the generalised Meyer context, with scaling functions with compactly 
supported Fourier transforms, allows the appearance of nonseparable 
wavelets, the context is nevertheless a constraining one. Also, we have 
shown that the intertwining relation is an exactitude that rules out 
compactly supported higher rank scaling functions. It is plausible that 
the setting of frames and higher rank GMRAs, which relaxes this con- 
dition in some way, may allow for the construction of compactly sup- 
ported higher rank frames. Moreover it will also be of interest to de- 
velop this theory further for dilation representations of Z d on IR n other 
than the basic ones we consider. 

Even in theoretical articles such as this one it is customary to make 
a few remarks concerning the potential efficiency of any new species of 
wavelets or frames. In the current climate of diverse and burgeoning 
applications this hardly seems necessary. However, we remark that 
nature often presents pairs of partially independent features with their 
own scaling aspects. Space-time scales are one obvious source of this. A 
purely spatial context can be found in the statistical theory of textures. 
This theory is suitable for localised analysis through wavelets (see [B] ) 
and, evidently, many natural textures have a rectangular emphasis. 

Our account is self-contained with complete proofs apart from a few 
basic standard lemmas. In particular in Section we reprove the con- 
struction and formula of a uniscaled dyadic wavelet associated with an 
MRA, and in Section 7 we give the construction scheme for classical 
Meyer wavelets. We refer the reader to the books of Wojtaszczyk |17| 
and Brattelli and Jorgensen [2] which are excellent sources for diverse 
wavelet theory. 
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2. Higher Rank Multi-resolution Analysis 

In this section we formulate the structure of a higher rank multires- 
olution analysis although we focus attention on the biscaled case of 
Definition 12.31 We obtain the frequency domain identification of a 
function f(x) in Vy in terms of a filter function rrif(£) and the Fourier 
transform </>(£). As preparation for the construction of wavelets from 
a BMRA we obtain a filter function matrix criterion in order that a 
given set {/i, . . . , f r } of unit vectors in L 2 {R d ) should generate an or- 
thonormal basis under the operations of translation and bidilation. 



2.1. Preliminaries. We take the Fourier transform / of an integrable 
function f(x) on IR d in the form 



f(0 



(2tt 



-i(t,x) 



f(x)dx, 



and define / for / in L 2 (R d ) by the usual unitary extension of the map 
/—>•/. The domain of Fourier transforms is referred to as the frequency 
domain and for emphasis is denoted M d . For a matrix A e GLd(M.) and 
a point x in the time domain M. d write Ax the product with x viewed as 
a column vector and define the unitary operator Da '■ L 2 (R d ) — > L 2 (M. d ) 
by 



(2) 



(D A f)(x) = \detA\ 1 ' 2 f(Ax). 



Recall that A is said to be an expansive matrix if all eigenvalues of A 
have absolute value greater than one. We say here that A is a dilation 
matrix if in addition ACZ d ) C Z d , that is, if A is a matrix of integers. 

A commuting dilation pair (A, B) is a pair of matrices in GL^i 



such that, after a suitable permutation of the standard basis of M. we 
have 



(3) 



.4 



A' 

Jw_ 



B 



h- q 
B' 



where A' and B' are dilation matrices in GL P (M) and GL q (M). Note 
that T = A n B m is a dilation matrix if and only if n > 1 and m > 1, 
while the matrices A* and 5 J are wea/c dilation matrices in the sense 
that their eigenvalues are not less than 1. In Sections 5, 6 and 7 we 
confine attention to wavelets in L 2 (M. d ) for a fundamental dilation pair 



.4 



'al 0" 




'I " 


/ 


B = 


131 



where a, (5 > 2 are integers. 
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2.2. Biscaled Multiresolution Analysis. We now define a biscaled 
analogue of a multiresolution analysis associated with a single dilation 
matrix. 

Definition 2.3. A biscaled multiresolution analysis, or BMRA, with 
respect to the commuting dilation pair (A, B) in GL^R) is pair (0, 03) 
where G L 2 (R. d ) and 03 = {Vy : G Z 2 } is a family of closed 

subspaces of L 2 (M. d ) such that 

(i) V i:j C Vj, m , for i<l, j <m, 

(ii) Dy^ ,o = ^, for alU,j, 

(iii) (1^)^ ^ = M. 

(iv) U( M)eZ ^M = ^(^), 

(v) PijPi m = Pi m Pij for all orthogonal projections P st '■ L 2 {R d ) — > 

(vi) {0(x — fc) : k G Z d } is an orthonormal basis for Vo,o- 

It is helpful to make the following distinctions. 

If we have merely conditions (i), (iii), (iv) then we say that 03 is a grid 
of subspaces. If, additionally, fl V^i = V m ,n with m = min{i,k}, 
n = min{j,l}, for all i,j,k,l, then we say that 03 is a lattice of sub- 
spaces. Finally, if, additionally, the stronger commuting condition (v) 
holds,then we say that 03 is a commuting lattice of subspaces. Thus a 
BMRA is a commuting lattice of type Z x Z which is generated by a 
single function and the dilation operators. 

We begin the analysis of BMRAs by identifying the various frequency 
domain descriptions of functions in Vij. 

Let / G Vij, T = A { B j and let t = detT. Then / o T _1 is in V ,o 
and so there is an expansion 

f(T- 1 x) = Cf(k)<f>(x - k) 

in L 2 (M, d ) for some square summable sequence {c/- (/c)} fcgZ d. This se- 
quence provides the filter coefficients of / for the dilation T. Define 
the 27rZ d -periodic function rrif(£) in L 2 (R d ) by 

™/(0 = ~£ c f (k)e-^ k l 

k&Z d 

This function is the filter function for / associated with T. Note that 
/ o T -1 (£) = f(T*£) and on the other hand 

k&L d 
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Thus f(T*£) = rrif(£)(j)(£) and so the filter function determines / in 
TVo t o by means of the formula 

(4) f(0=m f (T*~ 1 OHT*~ 1 0- 

We may view the filter function rrif(£) as being defined on 2nT d , in 

which case one readily sees that 



||"VlU 2 (27rTd) = 




In view of our applications and for notational simplicity we assume 
henceforth that A = A*, B = B* and hence T = T*. 

Note that m/(T _1 £) is periodic by translates from 27iT7* d and that 
we have the following converse. If <?(£) is 27rTZ d periodic and has 
square integrable restriction to 27rT a! then the function / in L 2 (R d ) 
with /(£) = lies in Vij. In particular we have the following 

criterion for membership of Vo,o which we state explicitly as it will play 
a role in the proof of the lattice property of a BMRA. 

Proposition 2.4. Let (</>, 9J) 5e a BMRA. Then a function f in L 2 (R d ) 
lies in Vo,o and only if /(£) = <?(£)0(O / or seme 2^2^ periodic func- 
tion g which has square summable restriction to 2nT d . 

Consider now the cosets of TZ d in the abelian group Z d , of which 
there are t in number, say E , . . . , E±-\. (See Proposition 5.5 in |17| 
for example.) Representative elements T , . . . , r f _i of the cosets are 
often known as digits and comprise a set of digits for T. 

For < % < t — 1 define the i th translate function for the filter 
as the function 

(5) m}(0=m / (e + 2 7 rT- 1 r 4 ). 

Note that from the norm correspondence above it follows that if / is a 
unit vector then the row vector formed by the t translates m l j(^) is a 
unit vector almost everywhere on 27rT d . 

Definition 2.5. Let (0, QJ) be a biscaled multiresolution analysis, let 
G Z, 2 , and let JF = {/ , . . . , f s } be an ordered set of functions in 
Vij. The (translation form) filter matrix for T is the function matrix 
Ut{£.) where 

where t = det A l B j . 

We remark that in Section 4 we define a related filter matrix L^r(£) 
whose rows are formed by the partial sums over cosets of the filter 
functions in the first column of Ujc. 
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The key filter matrix lemma below will be used in the construction 
of biscaled wavelets and in this connection it will be applied to various 
pairs X, Y from the spaces Vo,o, Vi t o, Vo,i, V\,\- The proof is similar to 
well known constructions in the the monoscaled case. 

Let X C Y be closed subspaces of L 2 (R d ) with Y = D T X and let 
be a scaling function for X, Y in the sense that the translates of 

form an orthonormal basis for X. In the light of our applications 
later we make the simplifying assumption that the digits are chosen 
so that the set T[0,27r] d is the essentially disjoint union of the sets 
[0,27r] d + 27rr i ,z = 0, ...,t- 1. 

First we note the following well known simple property of scaling 
functions which is equivalent to the orthonormality of translates. This 
will feature in the proof of Lemma 12.71 and in the construction of higher 
rank Meyer wavelets. 

Lemma 2.6. The set of Z d -translates of a function in L 2 (M. d ) forms 
an orthonormal set if and only if 

k&L d ^ ' 

Lemma 2.7. Let (0, X, T, Y, {Tj}) be as above and let T = {/o, . . . , f s } 
be a finite ordered set in Y with filter matrix Ujr{^). Then the system 
{fi(x — k)} ke x d ,i=o,...,s i> s an orthonormal set in Y if and only if for 
almost every £ the matrix Uj?(£) is a partial isometry with full range. 
Furthermore the system is an orthonormal basis for Y if and only if 
s = t — 1 and for almost every £ the matrix is unitary. 

Proof. Let /, g be two functions in T . We evaluate the inner product 

1 = (f(x — ki),g(x — /c 2 )) for hi, k 2 E Z d . Since f(x — k) has Fourier 
transform e~ l ^ ,k ^ /(£) the unitarity of the Fourier transform implies 

(6) / = / f(0W)e- i{tkl - k2) dC. 

Thus using (HI), the substitution r] = T -1 £, and the 2nZ d periodicity of 
the filter functions, we have 

(7) 1=1 m / (T- 1 e)m g (T- 1 OI0(r- 1 e)re- i<e ' fel - fc2> ^ 

(8) =t f m f (v)mM\kv) I V^'* 1 -* 2 ^ 

(9) =Y\t ( m f (ri)'^^\$(n-2TTk)\ 2 e- i{Tr >' kl - k2) dr]. 
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k&L d V ' 

and so 



(2tt) c 



11) / = —— \ m f { V )m g { V )e- i{TriM - k2) dr] 



[0,27r] d 



thus 



(12) 1 = 7^1 m f (T^)m g (T^)e^ kl - k2) dC 

Considering translates by 27iTj, and using periodicity in the exponential 
factor we conclude that the inner product {f(x — hi), g(x — k 2 )) is equal 
to 



t-i 



( 13 ) J^yI E m'fF-^mXT-H) e^"^. 

Now the sum function in the integrand is not merely 2irT'Z d periodic 
(in view of its terms) but is 2irZ d periodic by virtue of being a sum over 
all translates. Thus the integral is the (k% — A^th Fourier coefficient 
of the sum function. In particular, with / = g we deduce that the row 
vector function 

[m°(£) m}(0 ••■ mjr 1 ^)] 

is a unit vector almost everywhere. Now the lemma follows exactly as 
in the monoscaled case. ■ 

2.8. Separability and Higher Rank. We define a general higher 
rank multiresolution analysis for a commuting r-tuple (A\, . . . , A r ) in 
GLrf(R) to be a pair (0,23), where <p e L 2 (R d ) and QJ = {V*. : k G 
Z r } is a family of closed subspaces satisfying the r-fold version of the 
conditions (i),. . . ,(vi) of Definition 12.31 The r-tuple is assumed to take 
the form 

A\ = A l © I,i 2 © • • • © ld r i ■ ■ ■ i A r = 1^ © • • • © Idr-i © A' r , 

where d\ + c?2 + ■ ■ • + d r = d and A\ G GL^ (M) are dilation matrices 
(with integer entries). 

The simplest way to create such a multiresolution analysis is as the 
tensor product 

(01 © • • • © 4>r, 9Jl ® • • • ® 2J r ) 
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of MRAs (0i,2Ji), with 23 fc 




(*) . 



: % £ Z > , where 




Here the scaling function = 0i <g> • • • ® r is separable in the sense 
that there is an evident partition of the coordinates of L 2 (R d ) with 
<f)(x) = 4>i(x)4>2{x) ' ' ' 4>r{x), where 4>i(x) is a function of the coordinates 
in the ith partiton set. 

We say that a higher rank MRA (0, 03) of rank r is purely separable 
if it is unitarily equivalent to such an r-fold tensor product. For r > 3 
there can be partial forms of separability, which we do not discuss here, 
while for r = 2 we shall simply say that a BMRA (0, 23) is separable if 
it is unitarily equivalent to a 2-fold tensor product decomposition. 

If (<f), {Vuj\}) is a BMRA then selecting the diagonally labelled sub- 
spaces V iti gives rise to the MRA (0, {V iti }). We may thus observe from 
known facts for MRAs that there are redundancies in the conditions 
of Definition 12.31 Thus condition (iii) follows from (i), (ii), (iv), (vi). 
See Lemma 17.11 In general however, we have not found that the pres- 
ence of MRAs in BMRAs provides any shortcuts to the construction 
of higher rank wavelets. One can develop formulae for filter functions 
of functions in Vk and pursue the theory for rank greater than two but 
we do not do so here. 

3. Construction of Wavelets from BMRAs 

We now show how to construct wavelets and wavelet sets from a given 
BMRA, and we elucidate the interrelationship between filter functions. 

Let (0,93) be a BMRA in L 2 (R d ) for the dilation pair (A, B). We 
refine some notation for filter functions as follows. If / lies in V^o then 
we write (£) for the 27rZ <i -periodic filter function arising from 

T = A, and if / lies in Vo,i we write mf(£) for the filter function for 
T = B. Thus if / lies in Vi,o H Vb,i then we have, from (j4j), 



/(0 



m A AA l mA'\) 



and 




(14) 0(AO = m£(O0(O; 

(is) 0(bo = <(O<KO- 
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Put B£ for £ in (jl4l) and use (|T5l) to obtain 

(16) j>(AB0=m$(B0m%(0m- 



Reciprocally 




Since 23 is a commuting lattice the subspaces Vi }0 QV 0}0 and V^iBV^o 
are orthogonal and so we may define 



w 0<0 = v hi © ((Yi,o © Vo,o) © (v 0> i e Vb, ) © vb.0) • 

Moreover, let = £>* 1 Z^W r ,o for (i, j) G Z 2 , so that 



Thus these spaces are orthogonal. Moreover since the intersection of 
the spaces is the zero space it follows that 



and since the union of the V id is dense, L 2 (R d ) is the Hilbert space 
direct sum of all the W m>n . 

As in the monoscaled theory, if an explicit orthonormal basis 
{ifti, . . . ,ip s } can be constructed for the subspace W 0t0 then this is a 
wavelet set in the sense of the following definition. 

Definition 3.1. Let (A, B) be a commuting dilation pair in GL^(R). 
Then {ipi, . . . , ip s } is a wavelet set for (A, B) if 



is an orthonormal basis in L 2 (R d ). 

We now show how one can construct a wavelet set in Wo,o by means 
of a nested Gram-Schmidt orthogonalisation process and repeated ap- 
plications of both directions of the equivalences given in Lemma 12.71 

It is convenient to introduce the following terminology which antici- 
pates the construction process. As before we consider a BMRA (0, 23) 
for the dilation pair (A, B), where 23 satisfies the lattice condition. 
Also p = det A, q = det B. 




(m,n)<(i,j) 
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Definition 3.2. A wavelet family for the BMRA (0, 03) is a set 

T = {iPt ^t-i, V>f , • • • , Vf-i, • • • , ^s] 

where s = (p — i)(q — 1), and {ipf} (respectively {ipf}, respectively 
{ipk}) is an orthonormal set whose 1? translates form an orthonormal 
basis for V^o Vo,o (respectively Vo,i V 0t0 , respectively H^o.o)- 

Definition 3.3. A filter bank for the BMRA (0, 03) is a set of functions 
T in L 2 (R d ), 

F = {(p,tpi, ■ ■ -,^-1,^1, ■ ■ ■ ^f-i.V'i, • • • ,^s} , 

where the functions ipf E Vi t o Vb,o, £ Vb,i Vo,o, and ipk G H/o,o 
are such that 

(i) the associated (pg — 1) x (pq — 1) filter matrix U^^) for T = AI? 
is unitary almost everywhere, 

(ii) the p x p (respectively q x q) filter matrix £/ A (£) (respectively 
U B {^)) for the set {0, ipf-, . . . , V^-i} and T = A (respectively 

{0, i^i , . . . , V'^-i} an d T = 5) is unitary almost everywhere. 

If T is a filter bank as above then by (i) and Lemma I27T1 the functions 
of J- and their Z d translates provide an orthonormal basis for V\,\- Also, 
by (ii) and Lemma I2TT1 applied twice, the set {0, . . . , ^fx) (resp. 

, . . . ,V'^-i}) has translates forming an orthonormal basis for Vi.o 
(resp. Vo.i). It follows that the set {ipi, . . . ,ip s } has translates which 
form an orthonormal basis of PVo.o an d so this set is a wavelet set. 

Although Theorem 13.41 is stated as an existence theorem, the proof 
provides a recipe for construction which we shall carry out in the next 
section. 

Theorem 3.4. Let (0, 03) be a rank 2 MRA with respect to the com- 
muting dilation pair (A, B). Then there exists a wavelet set for (A, B). 

Proof. From the preceding discussion it suffices to construct a filter 
bank T = {0} U T. 

From and the dilation T = A construct the row vector valued 
function of £ given by the normalised row of translated functions for 
m^(£). This has the form 

[<(£) rr4{£ + 27^-^) • • -mfc + 2nA- 1 d p - 1 )} 

where d±, . . . , d p -± in IP x {0}, together with do = 0, give a set of 
digits for A. Precisely as in the monoscaled theory we may apply the 
Gram-Schmidt process to any full rank pxp completion of this row (by 
rows which are similarly translates of their first entry) to obtain a pxp 
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unitary matrix- valued function £/(£). We thus obtain U(£) = Ug{£) for 
a family Q = {0, ipf, . . . , V'^-i} where each ipf lies in V^o V^o- By 
Lemma \2.7\ these functions are orthonormal, with translates forming 
an orthonormal basis of V\$ Vo,o- 

In a similar way, using 0, the dilation T = B and a set of digits 
in {0} x 7L q for B, construct a unitary matrix Utc(l;) and orthonormal 
set K = {<p,ipi, . . . ,ipq_i} for which the functions ipf have translates 
forming an orthonormal basis for Vq,i © ^o,o- 

Consider now the union, 

Tab = {0, if>t ■ ■ ■ , , • • • , 

Since 23 is a lattice the ^ and the ^ are orthogonal. Moreover 
elements of Tab have orthonormal translates in V\ t i Vo,o- It thus 
follows from Lemma 12.71 again that for T = AB the (p + q — 1) x pq 
filter matrix for Tab and T, denoted Up AB (£), is a partial isometry 
almost everywhere. 

As before we may complete Up AB (£) to a unitary pq x matrix 
which is the filter matrix of a family Tab U {^i, • • • , V's}- This is the 
desired filter bank and, by Lemma 12.71 yet again, {-01, . . . ,^> s } is the 
desired wavelet set. ■ 



4. Latin Square Wavelets. 

4.1. Coset Filter Matrices. First we introduce a companion matrix 
U'-piC) for the filter matrix Ujr(£,) determined by functions {/ 0) • • • , ft-i} 
in Vi j and a fixed dilation unitary T, as given in Definition 12.51 This 
companion txt matrix uses cosets rather than translates and is unitary 
if and only if Uj^{^) is unitary. Furthermore, and it takes a particularly 
simple form in the case of the Latin Square wavelets. We were unable 
to find a reference for this equivalence and so give the detail here. 

Let T be a (possibly weak) dilation matrix, let t = detT and let 
E , . . . , Et-i be the cosets of TZ d . Let do, ... , d t -i be representative 
digits, with do = and E = TZ 2 . We are interested in the case 
T = A % Bi with spaces Vo,o and Vij. For a function / in V{j we have the 
coefficients Cf(k) for / and T as before, determining the filter function 

fcez d 

Recall that the matrix Up(£) is determined by its first column which 
consists of the filter functions m^(£),0 < I < t — 1. The rows are 
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formed by the translates (£), 1 < i < t — 1. Consider the coset sum 

= 7 E c A (j>^ 



4= Y] Cfi( d P + T k)e~ 



y/i 



1 e -j(?+27rT- 1 ,i I ,d p } 



where /x/; (£) is defined as the bracketed sum. Then 
m} iiF (0 = m i i,(e + 27rT- 1 d i ) 

= ^(e + 27rT-^)i= 
by the T~ X Z 2 periodicity of /ff, (£), where 

Note that D = {Dj,k) is a unitary valued matrix and 

t-l i-l 
p=0 p=0 

Thus 

where U'-p, which we call the coset filter matrix for T and T, is defined 
by 



4.2. Latin Square Wavelets. Let A = [g §] , B = [J§] be (weak) 
dilation matrices providing a commuting dilation pair (A, S) and uni- 
tary dilation operators Da,Db on L 2 (IR 2 ). Let be the characteristic 
function of the unit square [0, l] 2 in IR 2 . Then and (A, B) generate a 
BMRA (0, 33) in L 2 (1R 2 ). In fact V is simply the tensor product BMRA 
of two copies of the triadic Haar wavelet MRA. We have 

2 2 

0((A5)- 1 x) = ^^0(x-(z,j)). 

The distinctiveness of this scaling relation is that 4>((AB)~ 1 x) is simply 
a linear combination of translates of by a set of digits for T = AB. 
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This property, as we shall see, persists in the wavelets that we construct 
for (</>, 93) by filter matrix completion method of Theorem 13.41 We have 



4>(A 1 x) = (j){x) + 4>(xi - 1, x 2 ) + 4>{xi - 2, x 2 ) 



and so 



l + e~ i?1 +e~ i2 ^ 



Let (0, 0), (1, 0), (2, 0) be the natural set of digits for A. Then the 
three coset functions for m4(^) are simply the pure frequency functions 

-?=, ^e - ^ 1 ; Ts 6 " 42 ^ 1, ^° com pl e t e the row 

1 



v/3 



[1 



-i2£ 



to a unitary matrix valued function of coset functions we first complete 

to a 3 x 3 unitary scalar matrix. One such completion gives the desired 
completion 

- i i i - 

Vs Vs 

11-2 

v/6 Ve V6 
J_ L n 

-V2 V2 U - 

where M^(£) = diag(l, e'^e - * 2 ^ 1 ). A similar completion matrix U' B (£i, £2), 
with Mb(£) = diag(l, e~ l &e~ t2 &), is associated with 

m?(6,k) = J(l + e-*+e-^) 



^(£i,£ 2 ) 



and the digits (0, 0), (0, 1), (0, 2). 

Rows 2, 3 of the completions above provide functions ipf-, ijj^ ^2 
such that the five functions {0, -0f , 1P2} are part of a wavelet 

family in the sense of Definition 13.31 Moreover, for this orthonormal 
set the coset functions for digits for T = AB provide a partial isometry 



1 


1 


1 


1 


1 


1 


1 


1 


1 


:s 


3 




:i 


3 




3 


:s 


3 


\V2 








|^ 






-|V5 






1^ 



















IV2 


\V2 






£V5 


-|V2 








IV6 


-S^6 





±V6 


-iV6 





|^ 


-1^ 






where we have ordered the columns according to the digit order 
(0, 0), (0, 1), (0, 2), (1, 0), (1, 1), (1, 2), (2, 0), (2, 1), (2, 2). 
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We may now complete to a 9 x 9 unitary matrix, 



3 


3 


3 


3 


3 


3 


3 


3 


3 


Y'2 
6 


6 


v-2 
6 


Y'2 
(. 


y/2 
f, 


V2 

e 


3 


V2 
3 


%/2 
3 


y/6 

e 


vTi 
6 


V6 
6 


VE 

6 


VE 

6 


Ve 

6 











Vz 
e 


v'2 
(i 


V2 
3 


v% 

(i 


v" 

(i 


V2 
3 


v'2 
(i 


v'2 
(i 


yl 
3 


V£ 
e 


V6 
6 





Ve 
e 


V0 
6 





o 


Ve 

6 





VTo 

6 


\/io 

30 


2^10 
15 


2VT0 
15 


VTo 
30 


VTo 

6 


VTo 

30 


vTo 

15 


VTo 

30 





VTs 

30 


VTs 
30 


2VT5 
15 


2^15 
15 





2^15 
15 


\/l5 
10 


VTE 
30 


VTs 

15 


VTs 

6 


VTS 

10 





VTs 

15 


VTE 
15 


VTE 
15 


VTE 
10 


VTE 
30 


VTo 

10 





VTo 

10 





VTo 

10 


VTo 
10 


VTo 

10 


VTo 
10 


VTo 

5 



Just as the dilated scaling function <p(A~ 1 B~ 1 C,) is a linear combina- 
tion of digit translates of 0, so too are the wavelets, fa, fa, fa, ^4 which 
are determined by the last four rows of the completion. We can con- 
firm and understand the orthogonality of these wavelets by arranging 
the coefficients of -rj^fa, ~^, % I J 2, Tf^fai 77J5^4 as ^ n the diagram. One 



readily sees that the construction creates in this way a quadruple of 
latin squares which are pairwise orthogonal, that is, have vanishing in- 
ner products. Such constructs are natural to study in their own right, 
and indeed may be used to provide wavelets which, as here, are entirely 
natural variants of Haar wavelets. 

In summary, the arguments above have led to the following theorem 
where Xij( x J v) denotes the characteristic function of the set [0, l/3] 2 + 
(i/3,j/3). 

Theorem 4.3. Let tpi, . . . ,fa be the functions on M 2 given by 
X01 ~ 4X02 - 4xio - X11 + ^Xi2 - X20 



fa 
fa 
fa 
fa 



30 

s/To 
\/To 



(5x 



00 



. . , fa be the functions on 
4X02 - 4xio - X11 + 5xi2 - 
X01 + X02 ~ 4xio + 4xn + 4x20 - 3%2i - X22), 
5xoi + 3xo2 + 2xn - 2xi2 - 2^20 + 3%2i 



2X2i - X22), 



(2xoo 
(Xoo - 



X22J 



Then the set {fa, . 

{ 3 in+m)/2 fa (3 m x + k u 3 n y- 
is an orthonormal basis in L 2 



X02 + Xu ~ X12 ~ X20 ~ X21 + 2x22)- 
. . . ,ip s } is a bidyadic wavelet set. That is, the set 
k 2 ) : (m,n),(k u k 2 ) e Z 2 , 1 < i < 4} 
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-4 


5 


-1 




1 





-1 


-1 


-1 


2 




-1 


4 


CO 


5 


-4 


-1 







-4 


4 











3 


-2 


-1 




-1 


-1 


2 


-5 


2 


3 







1 


-1 


2 





to 




1 





-1 



Figure 1. Orthogonal Latin squares for the biscaled 
wavelets ipi, . . . ,1^4. 



5. Dyadic Biscaled Wavelets and Filter Formulae. 

In this section we derive the filter formula, as given in the introduc- 
tion, that corresponds to the commuting projection lattice property of 
a BMRA. 

We start by reproving the well-known fact that an MRA (</>, {Vj}) in 
L 2 (M) for the dyadic dilation matrix A = [2] has an essentially unique 
wavelet V'o- It is determined by the necessary and sufficient condition 
that 

M)=™*,(£/2)fe/2) 
where the filter m^ for ipo for the dilation T = A is given by 



where v(£) is an arbitrary 2-7r-periodic unimodular function in L°°(M.). 
To see this note that the row matrix function 

["ty(0 m M + 1T )\ 
determined by the scaling funciton <fi has a 2 x 2 unitary matrix com- 
pletion of the form 
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for every function v(£) as above. Furthermore, every unitary comple- 
tion in M 2 (L 2 (2irT)) necessarily has this form. It follows from Lemma 
12.71 that each such function ip has orthonormal translates forming a 
basis for the difference space V\ © V , and hence that if> is a wavelet. 
Conversely, if if) 1 is a wavelet, then by Lemma 12.71 {(f), if)'} necessarily 
has unitary filter matrix for T = A and so if)' necessarily is of the same 
form as i/jq. 

Consider now a dyadic biscaled wavelet if) by which we mean a 
wavelet for a BMRA (0,93) in L 2 (M) for the dilation matrices 



A = 



"2 


0" 




"l 


0" 





1 


, B = 





2 



obtained from Theorem 3.6. More precisely, since (detA — l)(det£> — 
1) = 1 the proof of Theorem 3.6 shows that there exists a wavelet set 
which is a singleton if> whose filter function m AB (^) for T = AB is 
associated with a filter bank {(f), if> A , if> B , -0}.This filter, m AB {^), arises 
from the step-wise unitary completion of the row matrix function 

H B (0 m£ B (£ + 7r(l,0)) m AB ^ + n(0,l)) m?{t + tt(1, 1))] . 

By our previous remarks we can make explicit how if) (x) and if> B (x) 
may be defined. Let us do this and recap the construction process for 
if). 

Define a function ip A by specifying its filter m A A (£) for T = A to 
have the form 

m£i(O = e-*mJ(£ + 7r(l,0)), 
and likewise define ip B via its filter function for T = B given by 

< B (O = e- l?2 mf(e + vr(0,l)). 

As we have already observed above, the 2 x 2-translate filter matrices 
for {0, ip A } and T = A, and for {0, ip B } and T = B, are unitary 
almost everywhere. It follows from Lemma 12.71 that ip A (resp. ip B ) 
has orthonormal translates spanning Vi,o Vo,o (resp. Vo,i Vb,o)- 
Since these difference spaces are orthogonal, by the commuting lattice 
property of a BMRA subspace grid, we obtain, via Lemma \2.7\ a 3 x 4 
partial isometry valued filter matrix function Uj?(£) for T = AB. It 
suffices to complete this to a 4 x 4 unitary- valued filter matrix in order 
to obtain an explicit filter function m AB which then determines the 
desired wavelet if). 

5.1. The Commuting Lattice Filter Relation. We now examine 
more directly the orthogonality of the rows of the 3x4 filter matrix 
for T = AB and {0, ifj B }. ' 
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For convenience we assume that the support of <fi contains [— 7r,7r] 2 . 
We are now regarding ip A as a function in V\ t \. This will necessarily 
have a filter function m A1 l (£) for T = AB satisfying 

i) A {ABi) 



We have, 



i> A {A{BO) 



m A A (BO<P(BO = m A A (BZ)m»(Om 



and so, almost everywhere on the support of <p, we have 

3-^(6 + 7r,26)mf(6, 6 



(0 



mjA(B0<(0 



By the support assumption and the 27rZ 2 periodicity of the filters the 
identitity holds almost everywhere. Similarly, 



e-* 2 mj(2£i,£2 + 7r)m£(£i,£2) 
To compactify notation we suppress £i and £2 and set, 



7T,0 



^1 = <(2£i,6), = mf(2(fc + tt),6) = ^, 



^1 = ™£(2£i,£ 2 + 7r), 



and so on, with 5^ , B*^ , . . . similarly associated with (£1, £2)- The 
3x4 filter matrix for {0, ■^ j4 ,^' B } is the product DU of the diagonal 
matrix D = diag{l, e - ^ 1 , e - ^ 2 } and the 3x4 matrix, 

" 4 0,0 r>0,0 att,0 tjtt,0 /1 0,0 R 0,7r att,0 r>ir,pi 

^l^" !,! A l,2 n l,l A l,2 n l,\ A \,2 D \,\ 



-T7T,0 



0.7T 



4 R Ul71 



-4°'°R ,r ' 7r 
^■1,2^1,1 



^ 0,0-^0,77 .^0-^0,77 ^.0,^-^0,0 ,7r,7r-p0,0 

/ 1 1,1- D 2,1 /i l,l- tJ 2,l ~~ ^1,1^2,1 ~~ ^2,1 

The first row of {7 can be written in the alternate form, 

r 4 0,0 R 0,0 4 7T,0 R 0,0 4 0,7T R 0,7T 4 7T,7T R 0,7rl 

l/M, 1^2,1 /1 1,1- D 2,1 /1 1,1- D 2,1 ^1,1 ^2,1 J 

in view of the intertwining relations 



rn 



iB ((i,( 2 ) = «J=«, etc. 



We note that the unitarity almost everywhere of the filter matrices 
for {(f>, tp A } and {(f),if> B } implies that almost everywhere 

|40,0|2 1 I An,0\2 1 Id0,0|2 , n0,iri2 i 
\ A l,l\ +1^1,1 I = X ' +l-°l,ll = L 

Thus, the inner product of rows 1 and 2 of U is 

^0,0 ^7r,0| R 0,0|2 a t,0 a 0,0 1 r tt,0 1 2 , a0,0 *tt,0 i r 0,tt i2 ^7r,0 ^0,01 rj7r,7r|2 
A \,2 A \,2 1-°!,! I — ^1,2 ^1,2 I I + /1 1,2 /1 1,2 I 25 !,! I — A \,2 A \,2\ D \,\ I 
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= A$A#(|Bj?l a + l^l 2 ) - A%A°$(\B$\* + \B%\*) = 0. 

Likewise, using the alternative form for row 1, the rows 1 and 3 of 
U are orthogonal. 

To recap, we have a 3 x 4 partial isometry translate filter matrix 
arising from the T = AB filters for (p, tp A , ip B , and 

m£?(0 = e-^m^(6 + vr,26K(6,6) = e^T^i 
mjffc) = e-*»mi(t u 6)^(2^6+^) = e^M?;^. 

As we have remarked, ^ A and ^ s lie in Vi t o Vb,o and Vb,i Vo,o 
respectively. This is a consequence of the commuting lattice property 
of a BMRA and it is for this reason that rows 3 and 4 are orthogonal 
almost everywhere. Conversely this necessary condition (together with 
the previous row orthogonality) is sufficient for the commuting lattice 
property. Indeed, the orthogonality of Vi,o Vo,o and Vb,i Vo,o follows 
from this (via Lemma 2.7 yet again) and the commuting projection 
lattice property follows. The row orthogonality is the following formula, 
which is written in expanded form in the introduction. 



ATXfi ,.0,0 R 0,0 R 0.7T ^0,0 A 1T,0 R 7T,0 R 0,7T 

I " /lO,7r R 0,7T R 0.0 , ,i0,0 47r,7r R 7r,7r R 0,0 r\ 

A l,2 A l,l n l,l n 2,l ' A l,2 A l,l n 2,l — U 



Remark 5.2. We have shown that for a dyadic BMRA scaling func- 
tion <fi there are two necessary conditions on the "marginal" filters 
m^(() and m?(£), namely the intertwining condition and the com- 
muting lattice (or orthogonality) condition above. In the final section 
we shall construct a function which satisfies these requirements and 
which defines a nonseparable BMRA. Evidently this nonseparable scal- 
ing function construction, ab initio, is considerably more complicated 
than that of constructing nonseparable wavelets, such as the LAtin 
square wavelets, from a given (possibly separable) BMRA. 

Remark 5.3. It would be interesting to pursue a "Riesz theory" of gen- 
eral not-necessarily-commuting lattices associated with scaling func- 
tions with, perhaps, Riesz basis translates. Going somewhat in this 
direction, we remark that the following is true (and the proof is rather 
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delicate). For a "noncommuting-BMRA", that is a pair (0, V) satis- 
fying all the axioms for a BMRA except the commuting lattice axiom 
(v), the subspace grid is necessarily a lattice. 

6. Impossibility of compact support for nonseparable 
bmra scaling functions 

In this section we show that if the scaling function of a BMRA (0, 03) 
in L 2 (M) is compactly supported then <f)(x,y) is separable. 

Lemma 6.1. Let a(£i, £2); &(£i> £2) be non-zero trigonometric polyno- 
mials with frequencies in 1? . Suppose that a, (3 > 2 are integers and 
for all £i,£ 2 , 

(19) a(£i,/36)&(6,&) = a(£i, 

T/ien a(£i,£ 2 ) = a(£i) anc ^ ^(^1^2) = ^(£2) / or some single variable 
trigonometric polynomials a(£i)> ^(£2)- 

Proof. Write 

Aff Afj 

(20) a(a,6)=£ $>j,*e iy€l+fc&) , 



and 



Aff Af| 



(21) &(&,6) = ££& 

where [L°, Mf ] x [L%, M 2 a ] =: Q a and [L\, M\] x [L*, M 2 b ] =: Q b are the 
minimal rectangles containing the support of the coefficients a^, bj^ 
respectively. Also define aj±, bj± to be zero outside their respective 
rectangles. For given p, q the (p, q)th term of a(£i, P&)b(£i, £2) is 



00 00 



(22) £ £ a p -j, k bj, q -i3k, 

j— 00 fc=— 00 

while the (p, g)th term of a(£i, £2)&(a£i, £2) * s 

00 00 

(23) ^ ^ ^ ^ &p—aj,kbj,q—k- 

j=—co k=—oo 

Consider the (p, /3M 2 a + Mf )th coefficient of a(6, /3£ 2 )&(£i, £2)- Since 



ap-j,kbj,q-f3k — a p-j,hbj,/3(M^-k)+M. 



1, 
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this term is nonzero only if k < M$ and (3(M% — k) < 0, that is, only 
if k = M%. Thus the Fourier coefficient is simply 

oo 

(24) a p _ jtI[ qb jtM b. 

j=-oo 

On the other hand, the (p,(3M% + Mf)th term of afa, 6)6(afi, f 2 ) 



is 

oo oo 



(25) a p-aj,kbj,(3M%+M%-k 



J 2 T 2 "2 

fc = — OO j = — oo 



There are nonzero terms in this sum only if k < Mf and flMg+M^—k < 
M\. Thus all terms, and the Fourier coefficient, are zero if Mf ^ 0. 

Assume, by way of contradiction that this is the case, so that, by 
the assumed identity a(£i, £2)&( a £i> £2) = a (£i ; (3£,2)b(£,i, £2) we have 

(26) ^ a p - jiM -b jtM b = 0. 

j=oo 

For the case p = Mf + Mf equation (126]) implies 

(27) a-M?,M$b M b M b = 0, 
hence at least one of Omj ,m£ , b M b M b is zero. Define 

(28) z a = max {j : a i Af « 7^ 0} , 

and 

(29) z b = max {j : &- M b ^ oj . 
Let p = z a + Zfe. We may now rewrite fT26l) as 

Zfr — l 00 

(30) ^ a p - jjM %b jiM b + a ZatM ?b z ^ M b + ^ a p-j,M^ ) j,Ml = °> 

J=-00 i = 2(, + l 

For j > z 6) = by the definition of ^. For j < z b , p — j > z a 

so by the definition of z a , a p -j t M% = 0. It then follows from (J30l) either 
a>z a) M% or ^.a/ 6 is zero, which is a contradiction, and so we must have 
Ml = 0. 

An analogous argument to the one just given, beginning with con- 
sideration of q = f3Li2 + L\, gives L 2 = and so 

(31) a(6,6)= E a ^ e ^' 
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Exchanging roles of the variables it follows that fe(£i, £2) is independent 
of £1, as required. ■ 

Theorem 6.2. Let (0,23) be a BMRA with respect to dilation pair 
(A, B) with scaling function £ L 2 (IR 2 ) of compact support. Then is 
separable. 

Proof. Recall that for a BMRA (0, QJ) with respect to dilation pair 
(A, B), for (£1,^2) we have the fundamental intertwining relation 

<(6,/36)™?(6,6) = <(6,6KKi,6)- 

If has compact support, then 0(£) is non- vanishing almost every- 
where. Furthermore 0(«6,6) and 0(6, /?6) are finite linear combina- 
tions of Z 2 -translates of and so the filters m^,m^ are trigonometric 
polynomials. By Lemma0m}fe, £ 2 ) = /(£i) and m|(£i,£ 2 ) = #(£ 2 ), 
where / and <? are trigonometric polynomials in one variable. It is 
routine to check that 0i(6) := 0(6,0), with filter /(£i) gives rise 
to a rank 1 univariate MRA with respect to dilation by a. Likewise 
02(6) '■= 0(0, £2) with filter #(£2) gives rise to an MRA with respect 
to dilation by (5. 

Using the filter relation (fl4l) N + 1 times gives 

= (ll/( a " n ^)) 0(«" n 6,6). 

\n=0 / 

As is compactly supported, is continuous. Furthermore / is a 
trigonometric polynomial, hence Lipschitz, and so the product 
n^=o/( a_n 6) converges almost uniformly, to -F(£i), say. Hence 

(32) 0(6,6) = F(6)0(O,£ 2 ). 
Similarly 

(33) 0(6,6) =G(£ 2 )0(£!,O), 

where (^(£2) = h m Af-»oo n^=o fi , (^~ n 6)- Hence, almost everywhere, we 
have F(£0 = 0(£ 1; O), G(£ 2 ) = 0(0, £ 2 ) and so 0(£ 1; £ 2 ) = 0i (6)0(6) 
almost everywhere as required. ■ 

7. Higher Rank Meyer BMRAs 

In this section we construct a family of bidyadic BMRAs which in- 
clude purely nonseparable examples. The construction is a higher rank 
version of the well-known method used by Meyer to construct wavelets 
belonging to the Schwartz class. In particular the Fourier transform of 
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the scaling function and the resulting wavelet have compact support. 
In fact the separable BMRAs obtained from the tensor product of two 
rank-1 dyadic Meyer type MRAs are included here as a special case. 
For our construction the scaling function and wavelet have discontinu- 
ous Fourier transforms which, as we have seen in the last section, are 
unavoidable. Thus our wavelets do not lie in the Schwartz class and it 
is not immediately obvious to what extent the decay may be improved. 

It is natural, by way of motivation and orientation, to recall the 
construction of Meyer wavelets, which we now do. 

Suppose that € L 2 (R) is a unit vector with orthonormal translates 
which satisfies the scaling relation 

(j){x/2) = ^a k (f)(x - k), 

with convergence in L 2 (1R), and is such that 0(£) is continuous at 0, 
with 0(0) 7^ 0. Then from and the scaling unitary for A = [2] one 
obtains an MRA (0, QJ). This fact is well-known ( [17] , Theorem 2.13). 

A scaling function of this type may be constructed by specifying its 
Fourier transform 9(0 = 0(£) to have the following three properties: 

(i) J2icz + 27rZ)| 2 = j-, almost everywhere. 

This condition is equivalent to the orthonormality of translates. (See 
Lemma 12.61) 

(ii) #(2£) = i/)(£)9(£), for some 27r-periodic function "0(C)- This is 
equivalent to the scaling relation above. 

(iii) 9(0 is continuous at with 9(0) ^ 0. 

To construct such a function 9(0 one may take the following route 
of Meyer and construct first a nonnegative function 9 on M. which is 
symmetric on [— 2tt, 2tt], with 

^)2 + ^_ 2vr )) 2 = _L on [ 0)27 r], 

2,71 

0(0 = -!=, for |e| < % 

6(0 = 0, for \0 > ^, 

Thus (i) holds, with at most two nonzero summands for each £. Let 
/(£) be the 27r-periodic extension of \/2~k6(20 for £ £ [— 7r,7r]. Then 
it follows that the scaling relation (ii) holds. If in addition 9(0 is 
continuous at with 9(0) ^ then the construction is complete. 
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It is completely elementary to construct a function 9 on K with the 
properties above. The main point in the construction is that, firstly, 
since #(£) = for |£| < we have the scaling relation 

/(£)£(£) = v^(9(20-^= = 0(2£), 

V Z7T 



which holds in fact for the bigger range |£| < 7r since #(2£), and hence 
/(£), are zero i n the range ^ < |£| < 7r. Thus, there is no obstacle 
to periodically extending /(£) to a function on IR and maintaining the 
scaling relation (ii). 

We are going to follow a similar procedure to construct a bidyadic 
scaling function in L 2 (M) which determines a multiresolution for A = 
[ o i ] 5 B = [qj]. However, while one can readily construct a nonnega- 
tive function #(£,£2) on [— 27r, 27r] 2 with the properties 

»(!) = ^fa!E(-j,j) 
»(0 = for ? £ |-~] 2 

and 

fl(£) 2 + #(£ " 2tt(1, 0)) 2 + fl(£ - 2tt(0, l)) 2 + 0(£ - 2vr(l, l)) 2 



4vr 2 ' 

there is no guarantee that one can periodically extend the functions 
^(0) ^ S (£) defined on the support of (9 by 

Our first main task is to construct 9 with extra structure so that 
this will be possible. This step is necessary because as we have seen in 
the general theory, if and A, B provide a BMRA then, by virtue of 
the subspace inclusions, will have periodic filter functions m^(() and 
(£) satisfying the intertwining relation. In fact we are arguing here 
in the reverse direction. We construct periodic extensions i/j a (£,), ip B {C)- 
These will be the filters for <p , since 

4>(M) = 0(2£i,£ 2 ) = 1^(00(0 = ^(£)fe), 

#(B£) = 0(£i,2£ 2 ) = ^ B (£)0(£) = ^(£)0(£). 
The intertwining condition follows from the equalities 

^(6,26)^(6,6)- /9(26 ' 26)U9(fl ' 261 



9(6,26); V 9(6,5: 
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^(£l,£2)^(2£l,£ 2 )- 



fl(2£) = ^(26,£ 2 ) \ ( 8(2Z U 2&) 
" 0(0 W(£i,6)y W(2£i,6) 

However, such a condition does not yet guarantee the orthogonality 
structure of the commuting lattice property and we must construct 9 
with further structure to ensure this. We do this in Theorem 7.2; we 
consider a function <p e L 2 (IR 2 ) and identify sufficient conditions on its 
Fourier transform (f) that ensure that <p, A, B determine a BMRA. 

Note first the following lemma. This follows from the rank one case 
(which is a basic fact; see |10] . [17j ). since every BMRA (<f), 23) contains 
the MRA {Vi,i :ieZ}. 

Lemma 7.1. Let (j) e L 2 (R d ) be such that {4>(x - k) : k e Z d } is 
an orthonormal set in L 2 (M. d ) spanning the closed subspace V^o- Let 
V it j = D' l A D j B V Qfi . Then fl^z* V itj = {0}. If, moreover, 0(0) ^ 
and </>(£) is continuous at 0, then U(j j)ez 2 ^m' ^ s dense in L 2 (R d ). 

We introduce notation for four families of rectangles which lie in the 
big square (— ^, ^f) 2 . 
For k e {0, 1}, let 



(M>) = ((-ir2 i+fc ^(-ir2 (i ^ )+fe |)x 

(-lV2 j+{1 ~ k) - f— lV'2 (1 ^' )+(1 ~ fc) — 
v ) 3 >v j 3 



and 



These rectangles lie between the big square and the small square (— ^-) 2 . 
In particular io,o,o is the rectangle which is the image of the north-east 
square (^, under A" 1 and the rectangles Iij,k have similar deter- 
minations. Note also that J^o is the closure of the union of the disjoint 

rectangles A~ k ((^-, ^) 2 ), k — 1,2, 

Also for i, j = 0, 1 let 

Z-47T , , -2lX , A /-47T , , -2ix , , 

#U = ( -3- + i(27r), -3- + i(27r)J x ^— + j(2n), — + j(2tt 

( -2l\ -7T A /-2-7T . , -7T . , 

,j = ( + + ) x I ~1T +J ^ , ~ +3 ^' 

Thus L 0i o lies in the south west corner of the small square. 

In the following theorem is assumed to be supported on the big 
square. Condition (d) is a condition on the restriction of to the four 
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Figure 2. Border rectangles in the support of 

corner squares (translates of ^) 2 ). Conditions (e) and (f) show 

how is determined on the border rectangles by the values of on 
pairs of corner squares. Condition (g) is an additional condition on the 
restrictions to corner squares. 

Despite the detail in conditions (a)-(g) the construction of examples 
of such functions is quite elementary. (Also there is further flexibility 
to arrange to be real- valued.) Indeed one may define on the four 
corner squares to comply with (d) and (g), and then is constructed 
(and uniquely determined) by the conditions (a),(b),(f). For example, 
in Figure 3 we show the regions of constancy of a function which 
takes constant values on the trangular subsets of the corner squares. 
Condition (d) is elementary and (g) holds trivially with both products 
in (g) identically zero. It is evident, from the triangularity of support 
in the corners, that and hence are not separable. 
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Figure 3. Constancy regions of a simple <fi example. 

Theorem 7.2. Let cf> G L 2 (IR 2 ) satisfy the following properties : 

(a) For (6,6) e M 2 we have < 0(6,6) < ^- 

(b) For (6,6) e (-f , f ) 2 we W 0(6,6) = 

(c) For (6,6) ^ (-f , f ) 2 ^ e W fe>&) = 0. 

(d) For (6,6) ^ (y> t) 2 ' ^(0 ^ s nonzero an d 

|0(6-2«,6-27rj)| 2 = i ^. 

i,36{0,l} 

fej For (6,6) e /(ij-fc) and k e {0,1}, 



l- #i(6,6) 2 



2tt 7 0i(6,6) 2 + ^(6,6) 



where 



12 h 



^1(6,6) =0(2 (1 ~ fc) 6,2 fc 6 
2 (6, 6) = 0(2 (1 " fe) 6 - 2vr(-l)^, 2 fc 6 - 2tt(-1)>'(1 - k)) 
(f) For (6, 6) e and z, j G {0, 1}, 

0(6,6) = 0(2 (1 " i) 6,2 t 6). 
fojfbr (6,6)e(f,f) 2 , 
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0(6, 6)0(6 - 2vr, £ 2 - 2tt) 2 - 0(6 - 2vr, 6)0(6, 6 - 2tt) 2 = 0. 

Let Vo,o &e i/ie closed subspace spanned by {0(6 — 6 — ^2) : ^2 £ 
Z}, Jet V-j = D^D'gVofl and let 2J = {V hj : i, j G Z}. T/ien (0,23) zs a 
BMRA with respect to the dilation pair A, B. 

Proof. We show that the conditions of Definition 12.31 hold. Let V^j be 
as above. Then part (ii) of Definition 12.31 is automatically satisfied. By 
Lemma [7.11 parts (hi) and (iv) of Definition 12.31 will follow from (vi). 
Let 

5(^1,6)= 10(6-27*, 6 -27Tfc 2 )| 2 , 

(fc 1 ,fc 2 )ez 2 

then (iv) will follow from Lemma [2.61 if we show that 5(6; 6) = 4^2 
almost everywhere. It is immediate from (b) and (c) of the theorem 
that this is the case for (6,6) £ ("F", if) 2 + 27r(/c 1 , fc 2 ). Likewise from 
(d) it follows for (6, 6) e (f , f ) 2 + 2 7 r(fc 1 , fc 2 ). For (£1,6) G / (0 ,o,o) = 
(f , ^) x (4f, f ) we have, from (e), 

5(6,6) H0(6,6)l 2 + 10(6,6 -2tt)| 2 
= J 0(26, 6) 2 

4vr 2 0(2^,£ 2 ) 2 + 0(2£i,6-27r)2 
| 1 0(26, 6 -2tt) 2 

47r2 0(6,6) 2 + 0(6,6-27r)2 

1 

as required. It is straightforward to carry out the preceding calculation 
for general I(i,j,k), hj, k G {0, 1} and so we obtain 

5(6,6) = 4^ fOT all (6,6) e {I (iJ , fc) + 2Z 2 : k G {0, 1}}. 

Now let (6, 6) e (0, f ) x (^, ^). Then there exists r G N such that 
f < 2 r 6 < f , and 

5(6,6) = 0(2 r 6,6) 2 + 0(2 r 6,6 - 2vr) 2 = J^. 

Again, this calculation may be repeated to show that 5(6,6) = 4~2 
for all (6,6) e + 2ttZ 2 : i, j G Z}. 

Hence we have shown that 5(6,6) — 4^2 almost everywhere, and 
so (0, 9J) satisfies (iii),(iv) and (vi) of Definition 12.31 

Next we show the BMRA inclusion property (i). From the definition 
of Vij it suffices to show VL^o C Vq,o and Vq ; _i C Vq,o- We show the 



30 



S. OLPHERT AND S.C. POWER 



former. By Proposition 12.41 this is equivalent to the equation 

0(26,6) = </(&, 6)0(6, 6) 
for some 27rZ 2 periodic function g with square summable restriction 



to 2%T . We show such a function exists. The equality holds trivially 

-4-7T 47T 

3 ' 3 



regardless of the value taken by g for (6,6) ^ ( — I 21 , t) 2 > nence if # 



may be constructed to be 27rZ 2 periodic on this square then its periodic 
extension will satisfy the equation everywhere. 

For (6,£ 2 ) such that f < |6| < f , 0(26,6) = but 0(6,6) ± 
so that we must take #(6, £2) = 0; also #(6,6) = #(6 ~~ 27r, 6) — 
for (6,6) such that -f- < 6 < -jf> so periodicity is preserved. For 
(6,6) e (~f,f) 2 we have 0(6,6) = ^, so we set #(6,6) = 
0(26,6); as 

((-f,f) 2 + 2vr(A ;i ,fc 2 ))n(-f,f) 2 = 

if (fci, ^2) 7^ (0, 0), periodicity is maintained. 

For the remaining (6,6) e (^, f ) 2 , both 0(6,6), and 0(26,6) 
are nonzero and we may define, as in the rank 1 case described earlier, 

0(26,6) 
<?(6,6) = -~— — -• 
0(6,6) 

Then for (6,6) £ ^(o,o,o) we have 



#(6,6) = 0(26,6) 



'0(26,6) 2 + 0(26,6-2tt) 2 



0(26,6) 



0(26,6) 2 + 0(26,6-2tt)2 
0(26, 6 -2tt) 



(0(26, (6 - 2tt) + 2tt)) 2 + (0(26, 6 - 2tt)) 2 



(0(26,6-2tt)) 2 



0(26, 6 -2tt) 

= — = #(6, 6 - 27r) 

0(6,6 -27T) 

so that # is periodic on /(o,o,o) U ^(0,1,0) ■ Again repeating a variation 
of this calculation, or by appealing to symmetry, we obtain periodicity 
for points in all I(%,j,k)- 

It remains to consider (6, 6) £ J(o,o)UJ(o,i). Periodicity of g for such 
(6,6) follows from the recursive definition of 0(6,6) 011 this interval 
and the periodicity of g on Iij t k- Hence g is periodic and VLi,o C Vo,o- 
By symmetry it follows that Vo-i C Vo,o and the other inclusions in 
Definition 12.31 (i) are satisfied. 
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Finally we must show that the spaces Vij form a commuting lattice. 
It is sufficient to show that (Vb,i Vo,o) -L (Vi,o © Vo,o) • As we have 
discussed in Section 5.1, we require 

(34) e- i « 1 ^/(6,6)=0, 

for almost every (6; 6) £ where 

7141,42) - ^1,2^1,1^2,1^1,1 - AAl^lA.l 

, /i7r,0 R 0.7r R 0,0 /i0,7r . ,1 0.0 R 7r,7r R 0,0 ,1 7r,7r 
~T ^1,2 - £, 2,1^ 4 -1,1 ' A \,2 a \,\ J3 2,l /1 l,l ' 

The factors here arise from the filters , of 0. Note that we have 
determined these filters as ratios, for example, m4 is given explicitly 
by the function #(656) above. 

Our aim is to show that the unwieldy expression above is equivalent 
to condition (g) for a function satisfying conditions (a)-(f). To that 
end we first observe that if on some rectangle at least one term in each 
of the products present in / vanishes, then / vanishes identically in 
that rectangle. 

For (£i,£ 2 ) G {K id } i>jem} , A?" = = = 5 2 °'° = 0, hence 
at least one term in each of the products in /(^i,^) is zer ° an d so 
/(6,6) = for almost every (£1,6) e {#i,jkj=o,i- For (£1,6) such 

4 7T,0 _ ^0,0 _ n Qrlr | Qr 

3 ^ |S2| ^ 3 



that < |£ 2 | < -jf , ^12 = ^12 = 0' an d so again for almost all such 



(6, £2) we have /(6, £2) — 0. Likewise, for (£1, £2) such that -y < |£i| < 



3 j -^2,1 = B21 = an d so /(6, £2) = for such (£1, £ 2 ). Hence we may 
restrict our attention to (£1,42) £ (— ^f, ^f) 2 - For (6>6) i n the central 
square (_|, A «fl = Al'J = = fljf = = B°J = 0, and 
hence /(6,6) = 0. For (£i,£ 2 ) such that f < |f 2 | < f and |^| < §, 
= A\f = Al'J = 0. Likewise for with f < |^| < f and 

|£ 2 | < f j -Bi'i = B^'i = B®'J = 0. Hence / vanishes for almost every 
(6,6) outside L i;i6{0 ,i}. 
Before proceeding, we observe that, for almost every (6)6) £ 

-f/r" I t \ A 2tT,0 7~)7T,0 T~)7T,7r ,177,0 4 7T,0 r>27T,0 r)7r,7r l27T,0 

/(6 + ?r,6) =^i,2 £i,i£ 2 ,i - A^B^ B 2 ' A A lt { 

i ^i27r,0 nir,if R 71 "^ /I 71 "! 71 " i /t""^ R 27r,7r R 7r,0 < 27r,7r 
+ ^1,2 -°1,1 D 2,\ A \,\ A l,2 n l,l - D 2,l /1 l,l > 

_ ,j0,0 R 7T,0 R 0,7T ,47r,0 171,0 R 0,0 R 0,7T 4 0,0 

— ^1, 2^1,1 - t5 2,l ^1,1 — ^1,2 D \,X D %\ 

1 /i0,0 R 7r,7r R 0,0 /jT,7r . ,i7r,0 R 0,7r r 0,0 ^0,7r 
+ - D 2,l /1 l,l ' A l,2 n l,l - D 2,l /i l,l ) 

= -/(6,6), 

where we have used the 2/tZ 2 periodicity of the filters (implicitly in 
y4°' 2 = A°' 2 and -B^f = -B 2 1). Similar calculations show that /(6> 6 + 
vt) = -/(4-vr, 6) and 6) = /(6+t, 6+tt), that is e-(&+&>/(6, 6) 
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is a 7rZ 2 -periodic function, and so it suffices to check that equation ([34 
holds on L Q for it to hold on all of the L^j. 

For (a, 6) e -^o,cb we have, by the definition of As±\ 

(35) A?;? = ^1^ = 2^(26,6) 

(36) = ^" 2& > 



V(2^ 1 ,2&)' + ^.(2? 1 + 27r,2& 
At this point we introduce the notation 
(37) $^ = 0(2e 1 + a,26 + 6) 2 , 

for a, b = 0, 2vt so that 



$0,0 

(38) 



$0,0 _|_ $2tt,0 

The other scaled filters also simplify. For example 



1,2 <K6 + tt,26) V $ 27r >° 

— ^$0,0 _|_ $27r,27r_ 

In fact on Lo,o all the scaled filters in the formula for /(6> 6) are 
expressible in terms of evaluated at the four points (26,26) + 
(«27r, j2vt), z, j G {0, 1}. With the other filters similarly expressed this 
leads to the equality 



0,0 



- $ 



2tt,0 



($2tt,0 _|_ $271,271) ($0,2vr + $27r,27r) 
($0,0 _|_ $2tt,0) ($0,0 _|_ $0,2tt) 

'($0,0 + $0,2vr) ($0,2tt _|_ $27r,27r) 



($0,0 + $2tt,0) ($2tt,0 _j_ $2tt 



27T / ($ 27r ' + $ 2,r ' 27r ) ($ ' + $ 27r ' ) 

($0,2?r _|_ $27r,2 7 r') ($0,0 + $0,2tt) 



($0,0 _|_ $0,2tt) ($0,0 + $2tt,0) 



($0,2tt _|_ $27r,27r) ($27r,0 _|_ $27r,27r) ' 

We simplify the expression by taking out the factor 

#(6,6) := (a/($ ' + $°> 2,r ) ($°>° + $2^,0) ($0,2tt _|_ $27r,27T) ($27r,0 + $27r,27r' 
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By the definition of 0, the function g is nonvanishing almost every on 
L 0) o and so condition 034|) holds if and only if 

= $0,0 ( $ 2^,0 + $2^,2^ ( $ 0,27r + $»r,2^ _ $ 2^,0 ( $ 0,0 + $ 0,2^ ( $ 0,2*- + $2^,2^ 
_ $ 0,2^ ( $ 2tt,0 + $2^ ( $ 0,0 + $2^,0) + $ 2^, 2 7r ( $ 0,0 + $0,2^ ( $ 0,0 + $ 2^,0) 

which, after some routine algebra, simplifies to 

($°<°$ 2 ^ _ $ 2 ^o$°^) ($°>° + $ 2 ^° + + $ 2 ^) = o. 

We observe that 

($0,0 + $2^,0 + $ 0,27r + $2^,2^ = _1_ 

since on rewriting in terms of it coincides condition (d). Thus ( 1341) 
holds almost everywhere if and only if, for almost every (6,6) £ -^0,0, 

(0(6,6)0(6 - 2tt,6 - 2tt)) 2 - (0(6 - 2tt,6)0(6,6 - 2tt)) 2 = 0, 

as required. ■ 

In summary, we have the following theorem. 

Theorem 7.3. Let A = , B = [Jj], Then, with respect to dilation 
by (A,B), there exists a nonseparable real valued wavelet ifj, associated 
with a nonseparable BMRA (0, QJ) in L 2 (M. 2 ), for which ip has compact 
support. 

Example 7.4. We finish with some further examples. Let {Hj}j e {o,i} 
be a partition of {-?Qj}j,.je{o,i} sucn that 

(i) {Hj} ie { 0jl } is invariant under the group Q of translation by ele- 
ments of 27rZ 2 , when {K it j}ij e T 0) n is viewed as a subset of M 2 /47rZ 2 , 

(ii) if (6,6) e s , then (-6,-6) e Hi. 

Observe that if is chosen to satisfy conditions (d) and (g) of Theo- 
rem 6.2 on Hq, then by defining 0(£ 1; 6) = 0(~6, ~6) for (6, 6) e Si 
the function then satisfies (d) and (g) on the whole of K . Consider 
the simplest case with constant on each set K^j flH for i,j G {0, 1}. 

Thus set 0(6,6) = « for (6,6) e (i\o,o n H ), and 0(6,6) = <5 for 
(6,6) ^ (-^1,1 HH ) for some constants a, 5 > satisfying a 2 +5 2 < ^ . 

Define 0(6, 6) = P for (6, 6) e (#1,0 n H ), where 




and 0(6,6) = 7 for (6,6) e (^0,1 n H ), where 
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We have now denned on S and, as suggested above, define 0(£i, £2) = 
£1, —£2) for (£i,£2) G Si. Thus we have defined to satisfy (d), (a) 
and (g) of Theorem 7.2, and so, as before, <fr, and 0, are determined 
and yield a BMRA. Note that we have 0(£i,£2) = 0(— £1, —£2); and so 
the scaling function is real valued. 

Notice that takes at most 4 values on each set I(i,j,k) and these val- 
ues, together with the values a, /3, 7, 5 and j- are the only values taken 

by 0. In particular, is a finite linear combination of characteristic 
functions. 
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